# Load the ggplot2 library
library(ggplot2)

# Define the function for PSP
psp_function <- function(omega) {
  return(omega * (1 - 0.15) / ((1 - 0.15) * omega + 0.05 * (1 - omega)))
}

# Generate a sequence of omega values
omega <- seq(0, 1, length.out = 1000)

# Create a data frame with omega and PSP values
data <- data.frame(omega = omega, PSP = psp_function(omega))

# Create the plot in black and white, with axis ticks but no grid lines
p <- ggplot(data, aes(x = omega, y = PSP)) +
  geom_line(color = "black") +                          # Black line
  geom_hline(yintercept = 0.95, linetype = "dashed", color = "black") +  # Black dashed line
  labs(x = expression(omega), y = "PSP") +
  theme_minimal() +
  theme(
    panel.grid = element_blank(),    # Remove all grid lines
  )

print(p)

# Export to high-resolution PDF
ggsave("Output/Figure1.pdf", plot = p, width = 3, height